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ABSTRACT 
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Aims. We report new results for the massive evolved and main sequence members of the young galactic cluster 
Oh DBS2003 179. We determine the physical parameters and investigate the high-mass stellar content of the cluster, as 

well as of its close vicinity. 

Methods. Our analysis is based on ISAAC/VLT moderate-resolution (R~4000) infrared spectroscopy of the brightest 
cluster members. We derive stellar parameters for sixteen of the stellar members, using full non-LTE modeling of the 
52 1 obtained spectra. 

Results. The cluster contains three late WN or WR/LBV stars (Obj 4, Obj 15, and Obj 20:MDM 32) and at least 5 
Olf and 5 OV stars. According to the Hertzsprung-Russell diagram for DBS2003 179, the WR stars show masses above 
85 Mq, the Olf stars are between 40 and 80 Mq, and the main sequence O stars are >20 Mq. There are indications of 
binarity for Obj 4 and Obj 11, and Obj 3 shows a variable spectrum. The cluster is surrounded by a continuous protostar 
| formation region most probably triggered by DBS2003 179. 

. Conclusions. We confirm that DBS2003 179 is young massive cluster (2.51O 4 X ) very close to the Galactic center at 

fC) ■ the distance of 7.9±0.8 kpc. 

Key words. Galaxy: open clusters and associations, stars: Wolf-Rayet, stars: early-type, stars: winds, general-Infrared: 
^\ 1 general open clusters and associations: individual ([DBS2003] 179 

o: 

1. Introduction tified near the center of the Galaxy (e.g. Cotera et al. 1996, 

Figer et al. 1997, Figer et al. 1999), but how they connected 
Where and how are massive stars forming in the Galaxy? to norma i stellar c l uste rs forming in the disk of the Milky 
Until recently, this question had not instigated significant Way was not yet known . Moreover, while efficient for study- 
investigation. However, the past decade has seen a surge in ing 1-2 kpc from the Sun (Carpenter 2000), the 2MASS sur- 
studies using recently available infrared surveys to search vey did not probe the inner regioil) wne re star formation is 
for deeply embedded massive stars and massive clusters. In at its grea test, and would not detect massive clusters if they 
2001, at the time of the first data release of the 2MASS were there ( Ha nson et al. 2010). Nevertheless, astronomers 
survey (Skrutskie et al. 2006), numerous new cluster candi- have discovered very massive clusters, rivaling anything 
dates were being discovered (Dutra & Bica 2001; Dutra et known in the Loca i Group of ga l ax i es in our Milky Way. 
al. 2003; Froebrich et al. 2007) , but no one had expected to Several new infrared surveys, such as GLIMPSE (Benjamin 
find very high mass young clusters. A few seemingly pecu- 2 003) and VVV (Minniti et al. 2010, Saito et al. 2012), 
liar, but clearly very massive, young clusters had been iden- are p roviding dee per imaging and greater spatial resolu- 

tion, which has unleashed a venerable "cottage industry" 

Send offprint requests to: J. Borissova of independent groups working to determine the massive 

* Based on observations gathered with ISAAC.VLT and cluster content of our Galaxy (Homeier et al. 2003, Bik et 

VISTA of the ESO and the 8.2-meter GEMINI telescope within a i. 2005, Mercer et al. 2005, Froebrich et al. 2007, Davies 
observing programs 81.D-0471, 179.B-2002, and GS-2009A-Q- 
81. 
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et al. 2007, Messineo et al. 2009, Mauerhan, Van Dyk & 
Morris 2011, Borissova et al. 2011). 

These all show that the Milky Way has many high-mass, 
young clusters, as predicted by Hanson (2003) based on the 
massive cluster populations found in external galaxies with 
Milky Way properties. But why is finding, cataloging and 
determining the properties of these massive clusters so im- 
portant? These young massive clusters are critical to our 
understanding of massive star evolution. The upper region 
of the Hertzsprung- Russell diagram contains different types 
of emission-line stars, including Of supergiants, luminous 
blue variables (LBVs), 03If/WN6 and Ofpe/WN9 stars 
(also called hot and cool slash stars), B[e], and Wolf-Rayet 
(WR) stars. All these classes represent different phases in 
the evolution of very massive stars. However, although it 
is commonly accepted that very massive O stars ultimately 
evolve into WR stars, the details of the evolutionary con- 
nections in the intermediate stages are still poorly under- 
stood. Their formation and early evolution are hard to in- 
vestigate, mainly because this period of the star's life is 
very short. The most direct way to investigate them is to 
identify and study massive stars and especially WR stars 
within star clusters, since this enables us to study a coeval 
population with uniform metallicity. For this reason, the 
most massive members of young galactic clusters deserve 
special attention. 

The open cluster candidate [DBS2003] 179 (hereafter 
DBS2003 179) was discovered by Dutra et al. (2003) via 
visual inspection of 2 MASS images during their survey of 
embedded clusters in the areas of known radio and op- 
tical nebulae. The cluster is located in the direction to- 
ward the HII ionized region G347.6+0.2. Borissova et al. 
(2005; hereafter Paper I) obtained high-quality infrared 
(IR) photometry and confirmed that DBS2003 179 is a 
member of the large family of "newly" discovered heav- 
ily obscured clusters. Based on the photometric analysis 
and on the so-called "10 th brightest star method" (Dutra 
et al 2003), we determined some of the cluster parame- 
ters: Ay — 19, (m— Af)o=13.5mag, age of ~7Myr, and to- 
tal mass of ~5.5 10 3 Mq. It was the first indication that 
the cluster is young and massive. However, the photometry 
alone is not enough to derive an accurate distance to the 
star cluster, and the uncertainty of the tenth brightest star 
method can easily be as much as a(m— M)o~2.5 mag or a 
factor of ^10 in distance (see the discussion in Paper I). 
To determine more accurate physical parameters for the 
cluster, we obtained medium-resolution (i?=9000) iVband 
spectra with the Infrared Spectrometer And Array Camera 
(ISAAC) (Moorwood et al. 1998) on the ESO Very Large 
Telescope (VLT) at Paranal, Chile, of selected cluster mem- 
bers and reported seven O stars and an Ofpe/WN9 star 
(Borissova et al. 2008; hereafter Paper II). Based on this 
we have calculated a distance of D~7.9kpc and a cluster 
mass approaching 10 4 Mq. The age of the cluster was es- 
timated to be between 2 and 5Myr. Mauerhan, Van Dyk 
& Morris (2011) reported another WN8-9 star associated 
with the cluster (the identification number MDM32) as part 
of their search for new Galactic WRs using the method of 
infrared color selection, and cross-correlation with X-ray 
point source catalogs. 

They calculated a distance to DBS2003179 of 6.4 kpc, 
using the calibration photometry values for WN8-9 stars 
given in Crowther et al. (2006). In contrast, in the most 
recent catalog of Milky Way high-mass clusters, Davies et 



al. (2012), report a kinematic distance of 9.0 kpc to the 
DBS2003 179. Thus according to Davies et al. (2012) the 
cluster holds a rather unique location on the far side of our 
Milky Way and possibly aligned with the Galactic Bar. 

A more in-depth investigation, using new and improved 
infrared spectroscopy of key cluster members, would more 
firmly determine its mass, but also its distance and its pos- 
sible association with the Galactic Bar. In this paper we 
report some new results for the massive main sequence and 
evolved cluster members. 

2. Sample selection 

In Paper II we used medium-resolution iVband spectra of 
DBS2003 179 obtained with the ISAAC on VLT at ESO to 
classify seven stars showing emission lines. Four of them 
(Obj 1, Obj 2, Obj 3 and Obj 5) are classified as early 04- 
06 V type stars, Obj 6 is most probably a young stel- 
lar object (YSO), while Obj 8 and Obj 4 are 04If-05If 
and Ofpe/WN9, respectively. Since the cluster seems mas- 
sive, we searched for additional emission-line stars. We 
constructed the narrow-band Br7 vs. 2.34 /im diagnostic 
color-magnitude diagram shown in Fig. [T] The Br7 im- 
ages were taken with Persson's Auxiliary Nasmyth Infrared 
Camera (PANIC) at the 6.5-meter Baade telescope at the 
Las Campanas Observatory, while for the continuum we 
used a 2.34/im ISAAC acquisition image (for more details 
on the photometric observations and data reduction see 
Papers I and II). All magnitudes are in the instrumen- 
tal systems. The stars with the spectral classification from 
Paper II are labeled. The diagnostic diagram indeed shows 
large Br7-Cont excess for Obj 4 and Obj 8, the stars with 
spectroscopically confirmed Br7 emission. The strong Br7 
excess of Obj 8 is puzzling, because the object does not show 
such strong spectral emission. Obj 6 is too faint and is not 
visible in the continuum image. Using the color excess of 
Obj 8 as a cut-off limit, as well as the position of the known 
emission line stars on the diagnostic diagram, we selected 
six candidate Br7 emission line stars. They are given in 
Table [T]and are labeled as Obj 11, Obj 12, Obj 13, Obj 14, 
Obj 15, and Obj 20. Three additional stars, which fell into 
the long slit (120arcsec length) during the observations, 
namely Obj 17, 18, and 19, are also given (see Sect. 4). 
As pointed out, the Obj 20 has recently been reported by 
Mauerhan, Van Dyk & Morris (2011). 

3. Observations and data reduction 

Medium- resolution spectra of the DBS2003 179 were ob- 
tained at the 8.2-m Unit 1 telescope of the ESO VLT, lo- 
cated on Cerro Paranal in Atacama, Chile, with ISAAC 
(Moorwood et al. 1998). ISAAC employs a Rockwell Hawaii 
1024x1024 array. The pixel scale of the detector is 0.146 
arcsec pixel -1 . Our spectra were obtained with the slit 
width of Of! 6, giving a spectral resolution of approximately 
4000 (slightly more in the if-band, and slightly less in 
the if-band). The grating was centered on two positions, 
1.705 and 2.150 /jm in the H and iVbands respectively. 
Unfortunately, some technical problems made it impossi- 
ble to obtain the -ff-band spectra for one of the slit posi- 
tions (containing the stars Objsll, 12, 13, and 20). The 
broadband JHKs imaging of DBS2003 179 was also carried 
out with the Aladdin detector of ISAAC during the same 
nights. 
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Table 1. Positional and photometric data for spectroscopically observed stars in DBS2003179. 



Object 


RA(J2000) 


DEC(J2000) 


J 


Err J 


H 


Err if 


K s 


Err K s 




hh:mm:ss 


deg:mm:ss 


mag 




mag 




mag 




Obj4 


17:11:31.88 


-39:10:46.9 


12.08 


0.01 


10.14 


0.03 


9.18 


0.02 


Objll 


17:11:31.33 


-39:10:55.3 


13.40 


0.07 


11.46 


0.06 


10.47 


0.05 


Obj 12 


17:11:31.55 


-39:10:53.0 


14.49 


0.05 


12.53 


0.02 


11.63 


0.04 


Obj 13 


17:11:31.71 


-39:10:51.1 


13.69 


0.02 


11.65 


0.02 


10.72 


0.02 


Obj 14 


17:11:31.67 


-39:10:46.9 


13.72 


0.02 


11.73 


0.03 


10.76 


0.03 


Obj 15 


17:11:31.80 


-39:10:46.8 


12.50 


0.02 


10.53 


0.03 


9.57 


0.03 


Obj 17 


17:11:31.69 


-39:10:44.4 


14.90 


0.03 


12.85 


0.02 


11.93 


0.03 


Obj 18 


17:11:31.85 


-39:10:29.6 


15.46 


0.03 


14.00 


0.05 


13.18 


0.03 


Obj 19 


17:11:32.01 


-39:10:20.1 


14.69 


0.06 


13.02 


0.04 


12.20 


0.05 


Obj 20(MDM32) 


17:11:33.03 


-39:10:39.7 


12.37 


0.02 


10.43 


0.03 


9.19 


0.02 
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Fig. 1. Results from the narrow-band photometry of 
DBS2003179 at Br 7 and Cont(2.34itm). The labels indi- 
cate the stars with known spectral types from Paper II. 
The emission line star candidates are marked with red cir- 
cles and are labeled. 

The near-infrared, high-resolution (i?«50000) spectra 
were obtained over five nights in February-June, 2009 with 
the Phoenix echelle spectrometer at Gemini-South obser- 
vatory with service mode in the if-band spectral range. 
Phoenix uses an Aladdin 512 x 1024 InSb array as the detec- 
tor, and it is operable throughout the range of sensitivity of 
InSb, 1-5 microns. The grating is a 63.4 degree echelle with 
32 lines per millimeter giving a dispersion of about 2 10 5 per 
pixel. Since Phoenix covers a very narrow wavelength range 
(0.5% of the central wavelength), we chose three spectral 
settings centered on 2.104, 2.143, and 2.179 /im. 

For the spectroscopic reduction, we used a combination 
of the Image Reduction and Analysis Facility (IRAFj], 
ESO eclipse package (Devillard 2001), and the Gemini 
pipeline. The basic data reduction included the flat-fielding, 
sky subtraction, extraction of one-dimensional spectra, 
wavelength calibration, and telluric correction (for more 
details see Hanson et al. 2010). The images were reduced 
following the typical procedures for the infrared, which in- 
cluded flat-fielding, sky subtraction, alignment, and com- 
bination into a final image. The photometry was obtained 

IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation 




Fig. 2. Identification of the stars with near-IR spectroscopy 
listed in Table [T] This is an ISAAC Kg band image with 
field of view of approximately lxl arcmin. North is up, east 
to the left. 



using the DA0PH0T package within IRAF. The photomet- 
ric calibration was performed by comparing our instrumen- 
tal magnitudes with the 2MASS measurements. The errors 
were calculated, taking the photometric errors and errors 
from the transformation to the standard system into ac- 
count. These magnitudes and errors obtained for the spec- 
troscopic targets are given in Table [1] 

4. Analysis of the spectroscopic targets. 

4.1. Spectral classification. 

The spectra of candidate emission line stars given in Table [T] 
are plotted in Figs. [3] and 0] in K and H (where available) 
bands. The preliminary raw spectral classification was done 
using available catalogs of i-C-band spectra of objects with 
spectral types derived from optical studies (Morris et al. 
1996; Hanson et al. 1996; Figer et al. 1997; Hanson et al. 
2005), as well as from the spectral catalogs of Martins & 
Coelho (2007), Crowther et al. (2006), Liermann, Hamann 
& Oskinova (2009), Mauerhan, Van Dyk & Morris (2011), 
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and Davies et al. (2011). The most prominent lines in the K 
and ii-band region (Brackett7 (4-7) 2.1661 /im (Br7); Hen 
A2.188 (7-10); Hei A2.1127 (3p 3 P° - 4s 3 S, triplet); Hen 
A1.692 (12-7); Hei A1.702 (ls 4 d 3 D e - ls 3 p 3 P°); Hi A1.668 
(11-4) ; Hi A1.722 (10-4)) were compared with the template 
spectra from these papers. 





Obj 17 



Obj 12 



Obi 10 



Obj5 




IU I I ■ 



1.699 1.722 
Wavelength (/^m) 



2.10 2.12 2.1-1 2.16 2.18 2.20 
Wavelength ifim) 



Fig. 3. The ISAAC spectra of DBS2003 179 stars. The dif- 
ferent spectra have been arbitrarily shifted along the flux 
F\ axis for clarity. The dashed lines indicate the rest-frame 
wavelengths of the spectral features. The best-fitting mod- 
els (green) are overplotted on the observed spectra (red). 

The equivalent widths (EWs) of DBS2003 179 spectro- 
scopic targets are listed in Table [2 They were measured on 
the normalized continuum spectra, by the IRAF task splot 
using the deblending function. The errors of the EW were 
estimated considering the signal-to-noise ratio measured at 
2.134 /im (see Col. 2 of Table [5]), the peak over continuum 
ratio of the line (see Bik et al. 2005) and the error from the 
telluric star subtraction (which was estimated to be ~10- 
15% in the worst cases). The EW of the emission lines are 
negative, and the symbol "NP" means that the line is not 
present in the spectrum. 

We derived the radial velocities of all stars with the 
IRAF task fxcor, which uses a cross-correlation Fourier 
method. Since the stars show very different characteris- 
tics, several radial velocity standards were used: for ob- 
jects 1, 2, 3, 5, 12, 13, 17, and 19 we used HD15629, 
an OV5 spectral type star with V ra d=— 48kms~ 1 ; Obj 8 
was cross-correlated with HD15570, an 04If spectral type 
star with V ra d=— 15kms _1 ; Obj 6 and Obj 14 were mea- 
sured by comparison with DM+493718, a Be spectral 






1.699 1.722 
Wavelength (fim) 



2.10 2.12 2.14 2.16 2 18 2.20 
Wavelength (//m) 



Fig. 4. The same as for Fig. [3] but for the stars with strong 
emission lines. 



type star; all absorption line stars were compared with 
HD191639 (B1V, V rad =-7kms- 1 ) and HD167785 (B2V, 
V ra d=— 10.6kms _1 ). The spectra were taken from Hanson 
et al (1996, 2005) spectral atlases. It is well known that the 
WR stars show broad and blended lines, which complicate 
the calculations of accurate radial velocities, so to estimate 
their radial velocities we used the IRAF task rvidlines. 
The task measures radial velocities from spectra by de- 
termining the wavelength shift in spectral lines relative to 
specified rest wavelengths. The basic usage consists of iden- 
tifying one or more spectral lines, entering the rest wave- 
lengths, and computing the average wavelength shift con- 
verted to a radial velocity. Thus the user can select the 
lines, which means that for WR stars it is possible to reject 
the broad and emission lines, which are not suitable for RV 
measurements. The lines are centered using Gaussian fits 
and when available, the purest line in our K band spectra, 
Hen A2.188, is taken as reference. 

The membership of the stars is determined on the ba- 
sis of the radial velocity histogram (see Fig. [5]). Most of 
the stars are concentrated within the radial velocity inter- 
val (—100, —50) kms -1 . The distribution histogram is then 
fitted with a Gaussian function. The mean radial velocity 
as given by Gaussian approximation is — 74±16 kms -1 . 
The uncertainty is determined from Poisson statistics of 
the measurements and the error of the wavelength solution. 
Objects 18 and 19 have very different radial velocities, lie 
relatively far from the cluster center (Fig. [SI bottom panel), 
and most probably are back- or foreground stars. Object 20 
also has more than 3cr of the fit radial velocity and could 
be field or runway star (see the next paragraph). 
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Fig. 5. Top: The radial velocity distribution of 
DBS2003 179. The line represents the best Gaussian 
fit. Bottom: The radial velocities for our spectroscopic 
targets as a function of distance from the cluster center. 
The horizontal lines represent the 3er dispersion of the 
Gaussian fit. The error bars represent the random error in 
determining the RV for each star, the outliers are labeled. 



We defined the following preliminary spectral classifica- 
tion of our targets: 

Object 11 shows broad and complex emission of Br7 and 
Nin/Oiii+Hei (hereafter 2.116/im complex), while Hen 
2.188 /im is moderately weak in absorption (Fig. [5]) and 
its spectral classification will be discussed later in this pa- 
per. Objects 12, 13, 17, and 19 are quite similar, with the 
2.116/im complex in moderately weak emission and Br7 
and Hen in absorption. For Objects 12, 13, and 17 the Hen 
is stronger than Br7. These are typical characteristics of 
O-type stars not later than 08. Based on their visual K- 
band magnitude, we assigned luminosity class V to Obj 12, 
Obj 17, and Obj 19 and I/III for Obj 13. Object 18 shows 
only lines in absorption, and the Hen line is absent in K- 
band spectra. Our spectral classification is 08-09, lumi- 
nosity class V. Object 14 shows Br7 in weak emission, no 
2.116 /im complex, and Hen in if-band spectra, but in mod- 
erately weak absorption in iJ-band spectra. The object is 
bright and most probably a supergiant. Objects 4, 15, and 
20 (MDM32, Mauerhan, Van Dyk & Morris 2011) show 
strong emission of Br7 and the 2.116 /im complex, and the 
Hen is in emission in Obj 15 and Obj 20 and in absorption in 
Obj 4. Following Crowther et al. (2006) we classified Obj 4, 
Obj 15, and Obj 20 as Ofpe/WN9, WN9, and WN8-9, re- 
spectively. This method of spectral classification, especially 
in the near-IR is correct within two subtypes, so we adopted 
this error for our estimates. 



4.2. Quantitative modeling. 

One of the goals of this work is to better constrain physical 
parameters for the massive stars. To do so, we proceeded 
with quantitative modeling of the spectra presented here, as 
well as of the stars reported in Paper II, which are obtained 
with the same spectrograph, with similar set-ups and the 
same reduction. They were separated into two groups. The 
first containing stars with no evidence of emission at Br7, 
with the second containing the stars that we classified as 
WR and Of. 

As in Hanson et al. (2010), we used cmf gen code (Hillier 
& Miller 1998). Although the spectra do not show many 
features, fairly complex atomic models were used in order 
to treat the line blanketing correctly. In addition we used 
306 levels for Nin and 343 levels Oin, which allowed us 
to reproduce the Nin 2.115 and Oin 2.116 /im blend. The 
absence of spectral features in our spectral range does not 
allow any estimations of the chemical composition of any 
elements heavier than oxygen so we adopted solar abun- 
dances for all them. It is commonly assumed that the winds 
are clumpy. We ran models with different clumping factors 
but the computed spectra showed that the sensitivity of 
the spectral features to its change at the resolution of our 
spectra are below the noise. The modeling presented in this 
paper use the same treatment of the clumping as described 
by Hillier et al. (2003). Following these authors we adopted 
a value of f=0.1. For the rest of the parameters we adopted 
the following fitting procedure. 

1. Velocity law. The wind velocity was assumed to follow 
the usual "/3"-type law. For the objects with Br7 in 
absorption we used Voo=2850 kms^ 1 and j3 — 0.8, as 
in the grid used by Martins et al. (2005) and in our 
previous paper (Hanson et al. 2010). 

For the emission line stars we used Br7 profile to esti- 
mate the terminal velocity. We find that the emission 
line profiles are not very sensitive to the changes in (3, 
but the ratio of the intensity of Hi 1.68 and Hi 1.74 to 
Br7 are sensitive to it. Higher values of (3 increase Hi 
1.74/Br7. Unfortunately, only the spectra of Obj 4 and 
Obj 15 are suitable for such an analysis. We determine 
(3 = 2 for Obj 4. The spectrum of Obj 15 is compatible 
with (3 = 1. We also assumed (3 — 1 for the stars without 
H spectra or with no emission in Hi 1.68 and Hi 1.74. 

2. Temperature and gravity acceleration. The wind of each 
model was fitted to a hydrostatic atmosphere with 
log g and T e g, both defined at the radius at which the 
Rosseland optical depth tr oss — 2/3. The Hei 1.70 /im, 
2.112 and Hen 1.69 /im and 2.18 /im features were used 
as temperature diagnostics. The gravity acceleration, 
log g, was determined for the stars with absorption lines. 
In addition to the shape of the wings of the hydro- 
gen and helium lines, we used the intensity of the Nin 
2.115/Oin 2.116 /im blend. The latter is very sensitive 
to the changes in logg, which makes the determination 
of the chemical composition of N and O difficult (see 
next item). The strength of the Hen 2.18/iin line is also 
sensitive to logg. The line became stronger with lower 
logg, which is to be expected. Thus the temperature 
and logg cannot be determined independently, which 
increases the uncertainty of the temperature and logg. 
Based on the step of the grid and the correlation be- 
tween the parameters, we assumed that the tempera- 
ture is determined within 2000K, although some of the 
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Table 2. EW and radial velocities of the observed stars. 



\Jl.Jj Ut, L 


S /N 


HI11-4 


Hen 


Hci 


HI10-4 


Nm /Ottt-I-TTot 
in in/ vviiin^ J- J- 1 — i 


y i 


Hen 


KV 






l.OOO 


1 RQ9 


1 709 


1 799 


9 1 1 fi 
Z. ±10 


Z. TOO 


9 1 fi£ 
Z. loo 


km s 


Obj4 


197 


-3.2 


0.8 


-7.2 


-6.4 


-6.7 


-45.2 


1.1 


-75±13 


Obj 11 


87 










-4.9 


-10.5 


1.5 


-54±7 


Obj 12 


78 










-3.7 


1.2 


1.4 


-62±9 


Obj 13 


71 










-2.8 


0.7 


1.1 


-66±6 


Obj 14 


52 


0.9 


1.2 


1.1 


0.7 


-2.3 


-2.9 


0.6 




Obj 15 


99 


-1.2 


-0.9 


-1.6 


-2.2 


-8.0 


-26.6 


-2.2 


-96±17 


Obj 17 


68 


0.9 


0.6 


0.9 


NP 


-1.7 


0.5 


0.9 


-76±15 


Obj 18 


35 


NP 


NP 


1.9 


2.1 


NP 


3.7 


NP 


-239±6 


Obj 19 


38 


0.8 


NP 


1.9 


NP 


-2.4 


2.1 


1.1 


-31±8 


Obj 20 MDM32 


80 










-20.7 


-47.3 


-2.9 


-200±30 



stars allowed for better precision. The stars that were 
fitted with models from the grid have an error in logg 
within 0.25 dex. Objects 15 and 20 have spectra that are 
not sensitive to logg. We adopted values that give rea- 
sonable masses for these two stars. The obtained tem- 
peratures and gravity accelerations, together with their 
estimated errors, are given in Table 3, and logg is in 
cgs-units. 

3. Luminosity and reddening. 

We determined the luminosity from the JHK magni- 
tudes. First we obtained reasonable values for T e ff and 
log g from the fit of the normalized spectra. Then we 
converted the observed magnitudes to fluxes using the 
Spitzer magnitude-to-flux converter Q. The equivalent 
fluxes were calculated from the computed spectra using 
the 2MASS filters' transmission curves 0. We then cal- 
culated a scaling factor for the adopted luminosity and 
the reddening for each object from the least squares 
comparison of its observed and calculated fluxes. We 
assumed the Cardelli et al. (1989) reddening law with 
R=3.1 and fixed the distance to 7.9 kpc. The derived lu- 
minosities are given in Table 3. We estimated the error 
of 0.2 dex in log(L/L©), taking the errors in tempera- 
ture, the distance and the reddening into account. 

4. Mass loss rate. Br7 was used to determine the mass- 
loss rate. As for T e g and logg, the errors are difficult 
to determine, due to the dependence of the equivalent 
width of Br7 from several other parameters. We esti- 
mated that the error in A4/y/J for the emission line 
stars should be less than 20%. The stars with Brj in 
absorption were modeled with the adopted mass-loss 
rate. Once the luminosity of the object was determined, 
the values of M. used in its model were scaled using Eq. 
(4) from Morisset & Georgiev (2009) to the transformed 
radius (Schmutz et al., 1989), and therefore the strength 
of the computed spectral lines is conserved for the de- 
termined luminosity. The results are given in Table 3. 

5. Chemical composition. The ratio between the He and H 
lines in both the H and K spectra were used to fix Hc/H 
ratio. The Nm 2.115/Om 2.116 /an blend was used to 
estimate the nitrogen and oxygen abundances. As men- 
tioned above, the blend is more sensitive to changes in 
log g than to the changes in the chemical composition, so 
the abundances determined based only on K spectrum 
and only on that spectral feature should be taken as 
rough estimates. Our analysis shows that, for all absorp- 
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tion line stars, a solar oxygen and nitrogen abundances 
is acceptable. The Civ 2.07 jim line is present at the 
blue edge of our low-resolution spectra (Hanson et al. 
2010). Although this line could be used to estimate the 
carbon abundance, in most of the cases, it is in the be- 
ginning of the ISAAC spectral range, is thus very noisy, 
and does not allow any rigorous analysis. In addition, 
the line has the same high sensitivity to \ogg as the 
Nm 2.115/Oin 2.116 /im blend, and it is also sensitive 
to the temperature. In most of the absorption line stars, 
the observed Civ 2.07 is compatible with solar carbon 
abundance models except for Obj 5, for which the model 
predicts an observable line, but none is detected. 
The stronger winds of the emission line stars make their 
spectra less sensitive to the underlying photosphere and 
therefore to logg. The Nm 2.115 line is strong and dom- 
inates the Nm 2.115/Om 2.116 fim blend. We still kept 
the oxygen abundance for these stars to its solar values, 
but we determined the nitrogen abundance, which in all 
cases is increased by a factor of 5-10 with respect to the 
solar. The Civ 2.07 is detectable in the spectra of the 
Obj 3, Obj 4, and Obj 8. Our current model for Obj 3 
is for its "active" phase (see below) based on our new 
spectrum, which does not cover the Civ 2.07 region. 
The Civ 2.07 is detected in the spectrum of its "pas- 
sive" phase. If this star is Of?p (see the discussion), the 
carbon lines are expected to vary, so the current model 
cannot be applied to its old spectrum, and we cannot de- 
rive the carbon abundance. Both Obj 4 and Obj 8 show 
much weaker Civ 2.07 than the model with solar car- 
bon abundance. To reconcile the observed and model 
spectra, we reduced the carbon abundance of the mod- 
els. The depleted carbon abundance and the increased 
nitrogen abundance are more evidence of the advanced 
evolution of these stars. 

The results are shown in Table [3] and model spectra 
are compared with those observed in Figs. [3] and [4] After 
analyzing all the information presented above we made 
the following conclusions for the most massive stars in 
DBS2003179. 

Object 3 shows strong emission in Br7 and the 
2.115/2.116/im complex, while the Hen is in absorption. 
With T cff =30000 K, it does not look significantly differ- 
ent from many mid-0 extreme stars, with an extended 
atmosphere, and it has most probably early 05/6If spec- 
tral type. Since the star was observed twice (in our 2007 
;tcqydir2C0ShISAAC runs), we compared its i^-band spectra 
and they showed significant differences. For example, in our 
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Table 3. Model parameters of the stars with spectra. 



Object 


Spec, type 


Toff 
[kK] 






l°g 9 
cgs-units 


log(L/L Q ) 


-M/v 7 / 

M Q /yr 


kms^ 1 






He/H 


log [N/H] 


Obj 1 


Oo/dV 


35.0 


±2 


5 


o c n i n n c 


5.54±0.2 




ZOZO 





8 


0.1 




-0.1 


Ubj 2 


uo/ov 


on r 
OZ.D 




r 



o.zoztU.zo 


O.Z4±U.z 


l.U / hj-Uo 


ZOZO 


n 
U 


o 
O 


U. 1 




0.0 


Obj 3 


UO/Dll 


on a 


±1 


o r 
2!) 


o i a_i_a i o 


O.o2±0.2 


O KOT7 1 ac 


OAA_L_ C A 


1 


U 


a oc_i_n n9K 
U.2o±0.U2o 





7±0.1 


al: Ad 

Ubj 4 


Of*™ MA7TVTA 

utpc/ wiNy 


oU.U 


±1 


O K 

Zo 


Z.yUztU.lz 


i k_i_a 9 
O.lOzhU.z 


Q OA T7 1 AC 


C AA-L- cn 

OUUzhoU 


i 

1 


n 


U. 1 


1 


0±0.1 


Obj 5 


UO/Ol 


or- A 

oo.U 


±2 


r 

a 


O.ZOztU.ZO 


O.42±0.2 


1 A '7T7 1 AC^ 


nonrfl 
ZOZO 


n 
U 


o 
8 


0.1 




0.0 


nu : cd 

Obj 6 


Yb(J: 


on r 

oz.o 


i - 

±0 


n 


o cnj_n cn 


4.02±0.0 


O OOT7 1 AT 


QAA I C A 


i 

1 


U 


0.1 




0.0 


i od 
Obj 8 


U4/D1 


o a a 


1 1 

±1 





o i a_i_a i k 


O.O2±0.2 


O OTT7 1 AC 


1 O K A_l_ K A 


1 


u 


i n_i_n A9K 
1.0±0.U2o 


U 


fi-l-A 9 


Ubj 1U 


Uo V 


oo r 

oo.o 


±1 


r 



o nn_LA i 
O.UU±U.O 


a aq_la o 
4.98±0.2 


OCT7 1 AC& 


22 ((J 


U 


8 


A 1 
0.1 




-0.1 


\J DJ ±Z 


V 


oo.u 


4-9 
ztZ 


r. 



O.OUztU.zO 


K /1 7-1-0 9 
0.4 / ztU.Z 


l.OUl-j-UO 


Zo^O 


n 
u 


o 


n i 

U. 1 




-0.1 


Obj 13 d 


04III 


34.0 


±1 


25 


3.10±0.12 


5.76±0.2 


1.05E-05 


1250±50 


1 





1.0 





6±0.1 


Obj 14 d 


05/6H 


32.5 


±2 





3.50±0.30 


5.72±0.2 


1.18E-06 


300±50 


1 





0.1 




0.0 


Obj 15 d 


WN8-9h 


37.0 


±1 





3.60±0.10 


6.29±0.2 


4.05E-05 


1000±50 


1 





0.15±0.025 


1 


1±0.1 


Obj 17 c 


06/7V 


35.0 


±2 


5 


3.50±0.25 


5.35±0.2 


1.29E-05 6 


2825 a 





8 


0.1 




-0.1 


Obj 18 c 


08/9V 


35.0 


±2 


5 


3.50±0.25 


4.96±0.2 


2.18E-06 6 


2400 a 





8 


0.1 




-0.1 


Obj 19 c 


05/6V 


35.0 


±2 


5 


3.50±0.25 


5.26±0.2 


l.HE-05 b 


2825 a 





8 


0.1 




-0.1 


Obj 20 d 


WN8-9 


35.0 


±2 


5 


3.40±0.25 


6.23±0.2 


1.18E-04 


1250±50 


1 





1.0±0.025 


1 


0±0.2 



a Value adopted in the grid of models. 6 Value scaled from the adopted in the grid to conserve the transformed radius (see text). 
c Models selected from the grid. d Models calculated to fit the particular spectrum. 



2008 ISAAC run the Br7 line is in emission (Fig. [SJ upper 
panel); as opposed to the 2007 run, where it is in absorp- 
tion), some variations in the 2.116 /im complex can also be 
seen. If confirmed, it would be an interesting case of a star 
in transition. The star resembles HD 191612 (Walborn & 
Howarth 2007) and the class of Of?p stars (see Naze et 
al. 2010 and references therein), which show similar vari- 
ations in the optical bands. The Of?p phenomenon covers 
stars presenting spectral variations (in the Balmer, Hel, 
CIII, Silll lines), strong CIII (A 4650 A) in emission (at 
least in some phases), narrow P Cygni or emission compo- 
nents in the Balmer hydrogen lines and He I lines (in some 
phases) , and ultraviolet wind lines weaker than for the typ- 
ical Of supergiants. Object 3 has increased He/H=0.25 and 
N/H=5 times solar (by mass), which closely resembles the 
composition of HD108 (Martins et al. 2010). The two stars 
have similar luminosities (logL/L Q = 5.7 for HD108 and 
log L/Lq = 5.5 for Obj 3), and although Obj 3 is cooler, the 
temperatures are compatible within the errors. It also has 
a peculiar position in the reddening-free color-magnitude 
diagram (see next paragraph). Thus we conclude that star 
Obj 3 is a candidate for another member of the Of ?p class 
(see Walborn et al. 2010, for the recent review of the five 
known members in the Milky Way of this type) however, 
long-term monitoring is needed to confirm this hypothesis. 

Object 4, classified as Ofpe/WN9 has a model temper- 
ature of T off =30000K, but the He/H ratio is low (0.1), and 
it shows a very high luminosity. It was also observed twice. 
As can be seen in Fig. [HI the 2.116 /mi complex does not 
change significantly (EW2oo7=— 6-0 vs. EW2oos=— 6.7), but 
in 2008 the Br7 line is stronger with EW2oo8=~ 45.2 vs. 
EW2oo7=— 37.2. This variability, if confirmed, could indi- 
cate presence of a close binary companion. Object 4, how- 
ever is placed very close to Obj 15, which was also classified 
as a WN8 star, and with 0.6 wide slit Obj 15 could have 
contaminated the spectrum. To check this, we searched 
for additional diagnostic lines and obtained high-resolution 
(R~50000) Phoenix, Gemini spectra in the two spectral 









_ CIV2.078 2.116 complex 


Brg 





2.08 2.10 2.12 2.14 2.16 2.18 2.20 

Wavelength (h/tj) 




2.08 2.10 2.12 2.14 2.16 2.18 2.20 

Wavelength (urn) 



Fig. 6. Upper panel: ISAAC, VLT, ESO, if -band spectra 
of Obj 3 taken in 2007 (thick black line) and 2008 (thin red 
line). Lower panel: The same but for Obj 4. 



settings centered on 2.104 and 2.179 /jm, respectively. In 
general, these spectral settings contain the following lines 
for the hot stars: Hel 2.1120 fim, CIII/NIII (8-7) 2.1032 ^m, 
HeII(23-9) 2.1786 ^m, Hel 2.1809 /Ltm, and Hel 2.1814 //m. 
In the upper panel of Fig. [7] are plotted three spectra that 
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CIII/NIII (8-7) 



2.100 2.102 2.104 2.106 2.108 2.110 
Wavelength (micron) 
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0.85 
0.80 



k 


'He2(23-9) 'Hel| He | 















2.176 2.178 2.180 2.182 2.184 
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Fig. 7. Phoenix spectra of Object 4 of DBS2003179. 

Table 4. Relative radial velocities of CIII/NIII line in 
Object 4. 



JD 


Central Wavelength 


RV 


CRV 


Day 


fim 


kins" 1 




54978.61149 


2.102911 





4 


54995.77851 


2.103142 


-42 


4 


55006.82398 


2.103249 


-62 


5 



were taken on 27 May, 13 June, and 24 June 2009, and 
centered on 2.104 /im. The lower panel shows the 2.179 /im 
setting taken on 4 June and 14 June 2009. The spectra are 
corrected with the mean radial velocity of the cluster. As 
can be seen from the figure, the CIII/NIII (8-7) lines at 
2.1038 /im are well defined and clearly present in emission. 
Moreover, the center and the shape of the line seem to be 
variable. The Julian date of the observation, center of the 
line (as measured from Gaussian approximation of the line 
profile by splot routine in IRAF), and the relative radial 
velocity with respect to the first spectrum and its error 
(measured by f xcor in IRAF) are given in Table 0] The 
width of the slit is small (0.34 arcsec) and is specially cen- 
tered to avoid the influence of the nearby Object 15. Thus, 
most probably Object 4 is a binary system. Assuming this, 
the relatively small variability of Br7 and the good fit of 
the spectrum with a single star model point to a relatively 
small contribution by the secondary to the observed spec- 
trum. Therefore one can conclude that the parameters de- 
rived from the spectra are representative of the brighter 
star in the system, but more data are needed to make more 
definitive conclusions. 

Object 5 was classified in Paper II as O-type dwarf star. 
Considering, however that the physical parameters derived 
for this star from the models are more compatible with the 
supergiants (even though under luminous) and its peculiar 



position on the color-magnitude diagram, we reclassify the 
star as 05/61. 

Object 6 has a K spectrum with only Br7 in emission, 
but the line shows two components. One is broad, corre- 
sponding to a stellar wind, and another much narrower 
which might originate in a circumstellar shell or another 
ionized gas on the line of sight. We tentatively classified 
this object as a YSO, but based only on one line, it is only 
a rough estimate. In the same way, we derived a tempera- 
ture for the star, such that the computed spectrum has Hei 
2.112 and Hen 2.18 lines weaker than the observed noise, 
and the mass loss rate and Voo fits Br7 profile, but the val- 
ues presented in Table 3 are only estimates. We reflected 
that by assigning larger error bars to its parameters. 

Objects 8 and 13 have a very uncommon spectrum. The 
spectra show a strong Hei/Niii/Oni complex with weak 
Hen 2.18 absorption line and very weak Bi'7. We repro- 
duced this spectrum with a helium and nitrogen-rich wind 
with depleted carbon. This makes the stars evolved, but 
their mass-loss rates are very low, which makes the objects 
unusual. Unfortunately, we only have if -band spectra for 
these objects, and the obtained set of parameters might not 
be unique. Other solutions with different chemical compo- 
sition might be possible. In any case, the stars have in- 
creased nitrogen abundance and are evolved from the main 
sequence. 

Object 11 has a very peculiar spectrum (Fig. [5J upper 
panel). The spectrum shows a double peak Br7, a P Cyg like 
Hei 2.11, and absorption in Hen 2.18. It is tempting to in- 
terpret the symmetric Br7 profile as the result of rapid rota- 
tion because the split between the two peaks is 540 kms -1 , 
equivalent to 270 kms -1 rotation. But a close inspection of 
the other two lines does not support this hypothesis. The 
lower panel of Fig. [8] shows the profiles of the three lines in 
the velocity space. It is obvious that Hei 2.11 profile is not 
P Cyg, but rather a composition of an emission line and 
an absorption component, centered on the same velocity as 
the Br7 central feature. On the other hand, Hen 2.18 is 
redshifted by ~40 kms -1 , and its profile is narrower than 
what one might expect for a 270 kms -1 rotational velocity 
(Fig. [51 lower panel). Based on that, one can suggest that 
the star is a binary with one component similar to Obj 4 
(Br7 and Hei 2.11 in emission and Hen 2.18 in absorption) 
and the other a cooler secondary with deep Br7 and Hei. 
We estimated L=4.34 10 5 L Q using Obj 4's model contin- 
uum as a reference and Ay=21.5. This makes the star not 
particularly bright, which is an argument against the bi- 
nary hypothesis. It seems that the star is peculiar and its 
interpretation requires spectra in wider wavelength range. 

Object 15 is classified as WN9, has a model temperature 
of T cf f =37000K, but the He/H ratio is low (0.15), and 
it shows a very high luminosity. This could indicate that 
the star is not a "classical" WR star, but rather a star in 
the transition 0->WN->LBV. Martins et al. (2008) found 
similar composition for the stars Bl, Fl, F6, and F9 of 
the Arches cluster. They classified these stars as WN8-9h 
even the composition resembles the one of a O star rather 
than the WR star. The authors conclude that the stars are 
descendants of the stars that are more massive than 60 M Q . 
The same conclusion can be applied to Obj 4 and Obj 15 of 
DBS2003179. 

Our model of Obj 20 is the most unsatisfactory. We were 
not able to reproduce the peculiar shape of Br7. The closest 
model has an increased He/H and N/H, high-mass loss rate, 
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Fig. 8. Upper panel: K spectrum of Object 11. Lower panel: 
from top to bottom Br7, Hel 2.11 and Hell 2.18 in velocity 
space. The lines are shifted by —70 kms -1 so the center of 
Br7 absorption component is at km/s. The two vertical 
dashed lines are at ± 270 km/s and mark the position of 
Br7 emission components. Notice that Hen 2.18 is narrower 
than Br7. 

and temperature adequate for a WN star. But the profile 
of Br7 has a peculiar rectangular shape which might be a 
result of a binary with two emission line components. The 
high luminosity, the temperature, and the chemical compo- 
sition of this star are similar to the parameters of HD5980 
(Georgiev et al. 2011). One can speculate that Obj 20 might 
manifest similar unstable behavior, and it needs further 
monitoring. 




J-K s 

Fig. 9. Top panels: Vector point diagrams of displacement 
of all stars within 35 arcsec of the cluster center (left) 
and spectroscopically confirmed cluster members (right) in 
units of ISAAC pixels (185 mas/pixel) after 5 years. Bottom 
panels: A zoom of the central 2x2 pixel vector point area. 
The most probable cluster members are within the red cir- 
cle with 0.25 pixels radius (left). They are marked with red 
crosses on the color-magnitude diagram of this zoomed area 
(right). All plots show only stars with root- mean-square po- 
sitions <0.08 pixels in each coordinate, the Objs 18, 19 and 
20 are also labeled. 



4.3. Relative proper motions. 

Objects 18, 19, and 20 are situated at the limit of the 
32 arcsec cluster radius determined in Paper II (see Fig. 0) 
with radial distances from the center of 20, 30, and 17 arc- 
sec, respectively, and thus they need special attention in 
order to verify their membership. We calculated the rela- 
tive proper motions of all stars within 35 arcsec of the clus- 
ter center. This small area was chosen in order to minimize 
the effects of distortion errors. We used the Anderson et 
al. (2006) method, which calculates the relative displace- 
ment between two epochs frames with respect to a local 
reference system of common stars. The best intrinsic astro- 
metric accuracy of 0.006-0.007 pixels in x and y coordinates 
measured by PSF photometry (the CENTER task) was ob- 
tained for a PANIC Br 7 image taken in 2003 and an ISAAC 
Xs-band image taken in 2008, which we chose as first- and 
second- epoch frames. At the distance of DBS2003 179 a 
five-year baseline is not enough to measure the proper mo- 
tion of the cluster stars from the ground. Instead, our intent 
was to track the possible high proper motion stars relative 
to our reference stars. The PPMXL catalog of positions and 
proper motions (Roeeser et al. 2010) lists 40 objects within 
a 1 arcmin radius centered on the cluster. We chose those 
25 of them with the smallest proper motions (<0.07 pixels 
in 5 years) to create our local reference system. The linear 
transformation between the two frames was computed using 
GEOMAP in IRAF, considering the shift, rotation, and scaling 



factor. Finally, we computed a global transformation-based 
displacement for each star defined as the difference between 
the Frame 1 position and the Frame 1 position implied from 
Frame 2 (based on the positions of the common stars) . The 
calculated errors were associated with centroid measure- 
ments for each star and the root-mean-square errors asso- 
ciated with the transformation. We then plotted the vector- 
point diagram (Fig. |H1 top- left panel) of the displacement of 
all stars within 35 arcsec of the cluster center, and the spec- 
troscopically confirmed cluster members (top-right panel). 
On the basis of spectroscopically confirmed members, and 
after taking their errors into account, we outlined a circle 
with a 0.25 pixel radius that isolated the cluster and field 
stars. The stars within this circle were considered the most 
probable cluster members (bottom- left panel). The stars 
Obj 18, Obj 19, and Obj 20 are notably outside this zone. 
In the color-magnitude diagram of the two zones shown in 
the bottom right-hand panel of Fig. [9] (open circles are field 
stars, red crosses represent cluster members while Objs 18, 
19, and 20 are blue asterisks) Obj 20 falls exactly in the 
WR zone of the cluster CMD (together with Objs 4 and 
15), while Obj 19 is at the limit and Obj 18 is more than 3 
<j from the mean cluster line. 

The spectroscopically calculated reddening of ^4y = 14.5 
for Obj 18 and 15.4 mag for Obj 19 are lower than the me- 
dian value (Ay— 18) of the rest of the spectroscopic tar- 
gets, while Obj 20 shows no differences. The last check we 
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Fig. 10. Upper panel: The reddening free color of 
(J-H)-1.77(H-K S ) vs. K S -0.78(J-H) magnitude. The 
solid line is a mean cluster line, while the dashed lines rep- 
resent 3cr deviation. Bottom panel: WISE 22^m image with 
density contours overplotted. The stars with spectra are la- 
beled. 



performed for these three stars was to note their position 
on the reddening-free color-magnitude diagram in order to 
estimate the effect of differential reddening. We used the 
reddening-free parameter Q =(J—H)—1.77(H—Ks), as de- 
fined by Negueruela et al. (2007) for OB stars. Figure [POl 
shows this reddening-free parameter vs. a corrected K$ 
magnitude. Object 20 has a peculiar position, which could 
be due to local dust overdensity and/or dust production. 
The first possibility is unlikely if taking the position of 
the star on the WISE 22^m image into account (Fig. ITU)) . 
Object 18 is far from the mean cluster line, while Obj 3 also 
has a rather peculiar position. 

After combining these results with the radial velocity 
measurements, we conclude that Obj 20 most probably is 



a cluster member, although its runaway status cannot be 
excluded, while Objs 18 and 19 are field stars. 

5. Cluster parameters 

Mauerhan, Van Dyk & Morris (2011) used the calibration 
photometry values for WN8-9 stars from Crowther et al. 
(2006), and assumed Obj 20 (MDM32) as a cluster mem- 
ber to calculate a distance to DBS2003 179 of 6.4 kpc with 
the associated error of 0.5-1 mag. In contrast, Davies et al. 
(2012) assumed that the cluster is physically associated to 
a nearby molecular cloud and its massive YSOs to derive a 
far-side kinematic distance of 9 kpc. Using the same tech- 
nique described in Paper II and the intrinsic color and ab- 
solute magnitudes from Martins & Plez (2006), we recal- 
culated the distance to the cluster. The average value of 
stars Objs. 1, 2, 3, 8, 10, 12, 13, 14, and 17 gives distance 
to the cluster (m— Af ) =14.5±0.3mag or 7.9±0.8kpc, in 
excellent agreement with the result obtained in Paper II. 
We did not consider the WR stars (Objs. 4, 15, and 20) or 
the stars with very peculiar behavior (Objs. 5, 6, and 11). 
The error is calculated from a standard deviation of the 
mean value, quadratically added errors from photometry 
and the variance due to an uncertain spectral classifica- 
tion. It is comparable within the errors with the kinemat- 
ical distance reported by Davies et al. (2012), especially 
considering that, in general, the kinematically determined 
distances are longer than those obtained from spectroscopic 
parallaxes. Thus the cluster is very close to the galactic cen- 
ter and could be on the other side of the galactic center. 

The K s vs. (J-K s ) and K s vs. (H-K s ) color- 
magnitude diagrams of DBS2003 179 are plotted in Fig.fTTl 
The potential cluster members (N=230) are selected by 
statistical decontamination as described in Paper I and 
the proper motion analysis described above. The the- 
oretical main sequence (Schmidt-Kalcr 1982) is shown 
and is corrected for the adopted E(B— V)=5A and 
(m— M)o=14.5±0.3mag. The stars with known spectral 
type are labeled, the zones of WR and OV class stars, re- 
spectively. In general, the Ks magnitudes of the objects are 
consistent with the luminosity class obtained from the spec- 
tral classification. Object 5 is slightly underluminous for its 
classification as a supergiant. 

The age can be estimated by fitting the observed 
CMD with solar-metallicity Geneva isochrones (Lejeune 
& Schaerer 2001) in combination with pre- main sequence 
(PMS) isochrones (Siess, Dufour & Forestini 2000). Starting 
with the isochrones set to the previously determined dis- 
tance modulus and reddening, we apply shifts in magni- 
tude and color until the fitting statistics reach a minimum 
value (i.e. difference in magnitude and color of the stars 
from the isochrone should be minimal) for both the main 
sequence and PMS isochrones. We get an "age envelope" of 
2-5 Myr. Figure[TT] shows the isochrone fits superimposed 
on the decontaminated CMD. As can be seen in the figure, 
however, for such a young star cluster, the main sequence 
isochrones are almost vertical lines in near-infrared, and it 
is hard to determine the precise age even when using the 
PMS isochrone set, especially in this case, when the pho- 
tometry is not deep enough. For this reason we determine 
the age using the HR diagram for the most luminous stars in 
the cluster by comparing the distribution of the stars with 
the latest isochrones and stellar tracks (with rotation) taken 
from Ekstrom et al. (2012). The temperature and luminosi- 
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J-K H-K 

s s 



Fig. 11. Left panel: The Kg vs. (J— Kg) color-magnitude 
diagram of DBS2003 179. All stars in the field of view are 
shown as light circles, while the cluster members are plot- 
ted with bold circles. The theoretical main sequence from 
Schmidt-Kaler (1982) is drawn with a solid line corrected 
for the adopted distance modulus of (m— M) =14.5 mag 
and E(B— V)=bA. Stars with known spectral types are 
labeled. The large circle and the triangle schematically 
mark the zones of WR and OV class stars, respectively. 
Right panel: K$ vs. (H—Ks) color-magnitude diagram. The 
solar-metallicity Geneva isochrones (solid lines) of 2 (green), 
3(red), and 5(blue) Myrs (Lejeune & Schaerer 2001) and the 
corresponding PMS isochrones (dashed lines) from Siess, 
Dufour, & Forestini (2000) are shown. 
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Fig. 12. Hertzsprung-Russell diagram for DBS2003 179 
(circles) and Arches (open diamonds, Martins et al. 2008) 
with evolutionary tracks from Ekstrom et al. (2012) over- 
plotted. The solid horizontal lines represent the mass tracks 
for stellar masses 120, 85, 60, 40, 32, 25, and 20A^ Q . The 
vertical solid lines represent 3 and 5 Myr stellar isochrones 
from the Geneva library. 



ties are taken directly from our quantitative modeling (see 
Table [3]) . The WR stars seem to fit the 3 Myr isochrone 
well, while the rest of the stars seem to group around 5 
Myr. A similar spread in age between O stars and WR 
stars was reported by Martins et al. (2008) for the Arches 
cluster and Liermann et al. (2012) for the Quintuplet clus- 
ter. They found ages of 2 to 3 Myr for the most luminous 
WN stars of Arches, while the O stars cover a range of 2 to 
4 Myr. They conclude that this might be because the most 
massive stars formed last in the cluster. Similar arguments 
are applied by Liermann et al. (2012) for the difference in 
age of the OB and WN stars in the Quintuplet. Thus, we 
adopt the age of the massive star population of 4±0.7 Myr 
for DBS2003 179. The stars, which show possible binary 
components, are not taken into consideration. The compar- 
ison with the isochrones from Girardi et al. (2002, "Padova 
isochrones" ) corresponds well with the derived from Geneva 
set age envelope. The error of this method can be estimated 
as ±2 Myr taking into account the photometric errors, er- 
rors of model fits, distance and reddening determinations. 

Luminosity functions (LFs) were created for the cluster 
in bins of 0.5 if-band magnitudes. To calculate the com- 
pleteness for each magnitude bin, we created artificial stars 
within the given magnitude ranges on the cluster image. 



The completeness was calculated by finding the recovery 
fraction of artificial stars per magnitude, using the same 
detection methodology we employed for our PSF photom- 
etry (see Sect. [3]). 

The LF was converted to the initial mass function (IMF) 
using the 4 Myr isochrone of Geneva assuming a solar 
metallicity. Since our photometry does not extend deep 
enough to detect pre main sequence stars, they were not 
a concern, and the stars within 5 a around isochrone were 
considered to be main sequence. The WR stars were also 
not considered in calculating the IMF's slope. 

Figure [H shows the IMF for DBS2003 179. The slope 
was calculated over the mass range limited by the verti- 
cal dashed lines. The upper mass limit was set to where 
the IMF slope starts to skew due to the effects of satu- 
rated stars, with the lower limit set to where the effects 
of incompleteness affects the slope. Completeness-corrected 
IMF is shown with dotted line. Uncertainties are also shown 
and were calculated from the small-number stochastic un- 
certainties from Gehrels (1986), where the upper limit is 
gi ven by 1 + V aN + 0.75 and the lower limit is given by 
\/<IN — 0.25. In our case, dN is simply the number of ob- 
jects within a given interval. The Gehrels equations give an 
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Fig. 13. IMF of DBS2003179. The slope is calculated be- 
tween the mass limits denoted by vertical dashed lines. 
The completeness-corrected values are shown as a dotted 
histogram. Overplotted as a red dashed line is the least- 
squares fitted slope (based on the completeness-corrected 
histogram) . 



approximation of higher and lower limits for each bin for 
low number statistics. Applying the equations to each bin 
gave the uncertainties. 

Our calculated slope of r=— 1.49+0.35, with uncertain- 
ties that propagate from those on the individual mass bins, 
is consistent with the Salpeter slope (r=— 1.35). This slope 
was used to extrapolate the masses of the observed main 
sequence stars to provide a total mass for the cluster. 
For extrapolation of masses below <O.5A^0 we used the 
Kroupa (1998) value of T=— 0.3. This gives a total mass of 
25300^7500^0 f° r ^ ne cluster. These large uncertainties 
on the mass are a consequence of the small mass-sampling 
region for the IMF slope. As a result the total mass that 
we present here should be treated with caution until deeper 
photometry can be analyzed. 

Finally, we show objects associated with the clus- 
ter star- forming complex. According to SIMBAD, within 
a field of radius=10 arcmin around DBS2003179, we 
can see: three HII regions, IRAS 17079-3905, GAL 
347.60+00.21, and GAL 347.6+00.2; four infrared sources, 
2MASS J17113303-3910400, MSX5C G347.5845+00.2123, 
MSX5C G347.6316+00.2131, and IRAS 17078-3910; five 
radio sources, GRS 347.60 +00.20, GPSR 347.616+0.152, 
GPSR 347.602+0.246, GPA 347.63+0.20, and GPSR 
347.629+0.149; ten masers; and three dark nebulae, which 
are indicators of ongoing star formation. They are overplot- 
ted on a VVV Kg band 15 x 10 arcmin image (Minniti al et. 
2010, Saito et al. 2012) in Fig. [T?] These sources are located 
on the periphery of the cavity surrounding DBS2003 179, 
clearly visible on GLIMPSE and WISE images. 

Although not cluster members, the Objs. 7, 9 (discussed 
in Paper II), 18, and 19 are OB stars, indicating a diffuse 
population of early main sequence stars around the cluster. 
To search for additional objects of the same type we used 
the Negueruela et al. (2011) approach. The reddening- free 
parameter Q =(J—H)—1.77(H—Ks) was used to separate 
the early-type stars from late- type main sequence stars and 
red clump giants, in four regions associated with the cluster 
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Fig. 14. Top panel: The VVV K$ image of 15x10 arcmin 
region around DBS2003 179 with known SIMBAD objects 
overplotted. Bottom panel: The WISE, GLIMPSE, and 
VVV true color images of 6x10 arcmin region. 



itself, and the star formation indicators (HII regions, in- 
frared sources, masers, and radio sources). Figure [TBI shows 
this reddening free parameter vs. a corrected K$ magni- 
tude for such selected areas, as well as a comparison field. 
Originally, Negueruela et al. (2011) defined Q <0.0-0.1 
as a separating value for early-type stars. In our case, we 
defined this value using the right partition of the spec- 
troscopically confirmed OB stars, namely Q <0.4— 0.5. As 
can be seen in Fig. [T5J there is no overdensity of these 
stars in Fields 2, 3, and 4 with respect to the comparison 
field. In contrast, Field 1 centered on DBS2003179 (with 
radius 1 arcmin) shows an overdensity. The spatial resolu- 
tion of VVV images did not allow us to measure most of 
the DBS2003 179 cluster members due to crowding. Thus, 
there may be a larger OB association around the cluster, 
most likely preformed at the same time and from the same 
molecular cloud as DBS2003 179. Interestingly, except for 
Obj 6, only two other YSO objects have been reported to 
date in the close vicinity (15 arcmin radius) of the cluster. 
Analyzing the color-color diagram of the VVV frame, we 
can see some stars with infrared excess and/or PMS stars 
(Fig. I16|). However, our photometry is not deep enough to 
investigate this in more detail. 

In summary, we can conclude that, according to VW 
images, there are no other younger and brighter than 
Kg— 18 mag star cluster and/or stellar groups in the close 
vicinity of DBS2003179. However, taking the total lumi- 
nosity of modeled O and WR stars (approx 1.6 10 7 L/L®) 
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Fig. 15. The reddening-free color of (J-H)-1.77(H-K S ) 
vs. Kg— 0.78( J— H) magnitude. The stars with spectra are 
marked with red filled dots, gray dots are all stars in the 
selected regions, the blue triangles represent the stars from 
the comparison filed. The solid line is a separation line be- 
tween early OB stars and late-type main sequence stars and 
red clump giants (see the text). 

into account, we should expect some continuous or burst 
protostar formation that is most probably triggered by 
DBS2003 179 and the surrounding OB association. Indeed, 
several dust emission regions are visible on WISE and 
GLIMPSE. They will be analyzed in the next paper in our 
series. 

6. Summary 

We report new results for the massive evolved and main se- 
quence members of the young galactic cluster DBS2003 179. 
The cluster contains three late WN or WR/LBV stars 
(Obj 4, Obj 15, and Obj 20 (MDM 32)) and at least five Olf 
and five OV stars. According to the Hertzsprung- Russell 
diagram for DBS2003 179 where the modeled parameters 
of all spectroscopic targets are plotted, the WR stars show 
masses above 85 Mq, the Olf stars are between 40 Mq and 
80 Mo, and the main sequence O stars are >20 Mq. There 
are indications of binarity for Obj 4, Obj 11, and Obj 20; 
Obj 3 shows a variable spectrum; and Obj 15 could be a star 
in transition. Using the newly added spectroscopic targets 
and photometry, the distance and mass of the cluster was 



Fig. 16. VVV observed (J—H) vs. (H—Kq) color-color dia- 
gram. Only stars with photometric errors <0.05 are plotted. 
The lines represent the sequence of zero-reddening stars of 
luminosity class I (Koornneef 1983) and V (Schmidt-Kaler 
1982). The reddening vector for Ay ~ 15 mag is overplot- 
ted, and the dotted lines are parallel to the standard red- 
dening vector. 



redetermined, giving 7.9±0.8kpc and 2.51O 4 .M , respec- 
tively. The WR stars fit 3 Myr isochrone well, while most 
of the O stars are placed around the 5 Myr isochrone. The 
slope of the IMF is consistent with the Salpeter slope. The 
cluster is surrounded by a continuous protostar formation 
region that was most probably triggered by DBS2003 179. 
There is also an indication of diffuse OB association and 
pre-main sequence stars in the close vicinity of the cluster. 
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